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ABSTRACT 


The conventional approach to ship steering is to regard the ship as a single input, 
single output system without cross-coupling or interaction between speed, yaw and roll. 
This approach has found successful application, particularly in conventional vessels 
where the amount of cross-coupling is normally slight. But, as a result of tight 
maneuvering, the modern warship suffers severe cross-coupling effects because of large 
control surfaces, high speed and low tonnage. Consequently, the adoption of a 
multivariable approach to ship steering would appear to be more suited for the design 
of a steering control system. 

This thesis describes the results of a simulation study of pre-compensator design 
to suppress the undesirable cross-coupling effects between speed, yaw and roll. 

Simulation studies using DSL and Function Minimization are the basis for 
accomplishing the design. 

Simulation results presented indicate that the adoption of a multi-input, multi- 
output approach would result in a significant improvement in the combined steering 


and stabilization problem of a warship. 
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I. INTRODUCTION 


Modern warships must be highly maneuverable to satisfy numerous operational 
requirements. However, the steering characteristics of a modern ship are nonlinear, so 
severe interaction or cross-coupling exists between the control surface inputs and the 
controlled outputs. Therefore, warship steering is a complex multivariable control 
problem. 

Interaction between roll, vaw and speed are pronounced in warships because of 
their length to beam ratio and relatively large control surfaces. However, warships 
should be able to execute high speed maneuvers while maintaining their fighting 
capability. This 1s generally not possible due to severe cross-coupling. 

In this thesis, the development of a pre-compensator to reduce the cross-coupling 
effects which are present in the steering characteristics of a modern ship is introduced. 
The design method adopted by Roberts [Ref. 1] uses the Direct Nyquist Array (DNA) 
frequency response technique as defined by Rosenbrock [Ref. 2] and Fricker (Ref. 3]. 
This thesis uses the basic ship model and the results which were reported bv Roberts 
Which have the potential to produce improved seakeeping and ship stability. 
Emphasis in this thesis is placed on optimizing the parameters of a pre-compensator 
using Function Minimization (F.M.) via a digital computer. 

Simulation studies employed Function Minimization techniques together with the 
Dynamic Simulation Language (DSL) package. 


ha 


If. BASIC MODEL 


The model! and data used in the study is that proposed by Roberts [Ref. I]. The 
Structure of the ship model is shown in Figure 2.1 and the elements of the transfer 
function matrix, G(S), which were derived using curve fitting techniques to measured 


step and frequency response data, are given in Table 1. 


RUDDER 





Figure 2.1 Multivariable Structure of a Warship Model. 


The nonlinear nature of the ship dynamics is demonstrated by the change in 


steady state gain parameters as shown in Table 2. 


12 


eae at 
ELEMENTS OF MATRIX FOR WARSIIIP MODEL [REF. 1] 


G,j(S) = a 
- 45° + 0.245 + i 
K,.(-8.57S + 1) 
Gi(5) = 3 2 
53.35? + 17.178° + 9.528 + 1 
G,,($) = ae 
ee $(12s’ + 32.25s¢ + 11.28 + 1) 
K,,(10S + 1) 
G,(8) = ———— 
(8) 240s? + 58S* + 26S + 1 
G,,(S) = ar 
Me 240S* + 34S + 1 
K 33 
G..($)"= 3 
oa 245 + 1 
G35) = G(5) = G3(8) = 0 


TABLE 2 
GAIN VARIATION [REF. 1] 





For ease of analysis the system inputs and outputs are used as defined below in 
Maple: 


TAREE 


MAX VALUES OF INPUT ANDIOUTLCED 





INPUTS MAX VALUE 


Yaw ee +120 


Figure 2.2 and Figure 2.3 show the time responses for the ship model when we 
use 10° demands for rudder and fin at 12 kts. 
As can be seen in Figure 2.2, there are pronounced cross-coupling effects in both 


roll and speed outputs. 
In this thesis, 10% demand for rudder at 12 kts will be alwavs used for 


convenience. 
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Figure 2.2 Step Response for Uncompensated Warship Model 


at Rudder Demand. 
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Figure 2.3 Step Response for Uncompensated Warship Model 


at Fin Demand. 
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Itt. PRE-COMPENSATOR DESIGN 


The aim of including a pre-compensator, K(S), is to decouple the interaction 
present in G(S), thus enabling the reduction of interaction or cross-coupling in the 
system. The action of the pre-compensator is to propagate the three input demands in 
such a way that each input affects its associated output only. 

This chapter was extracted from Roberts [Ref. 1]. 

The system configuration with the pre-compensator included is shown in Figure 
Dal. 


WANTED ROLL 


ANTED COURSE 


WANTED SPEED 





Figure 3.1 Compensated System Configuration. 


The approach used in this study is to develop a pre-compensator, K(S), which 
totally diagonalizes the pair G(S)K(S) ( = Q(S) ). This is in effect non-interacting 
control and can result in the elements of K(S) having high order which may prove a 
problem when it comes to implementation, particularly if this 1s to be achieved using 
analog techniques. If necessary the complexity of the elements of K(S) can be reduced 


using standard reduction routines while maintaining diagonal dominance in G(S) K(S). 
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The diagonalizing pre-compensator was produced using the method suggested by 
Fricker [Ref. 3]. This technique produces an initial “ideal” pre-compensator K(S) which 
diagonalizes Q(S). The elements of K(S), which can be high order, are then reduced to 
a simpler form while ensuring that diagonal dominance in Q(S) is maintained. This 


involves expressing G(S) as: 


G(S) = R(S) A(S) (ean 3.1) 


where R(S) is a diagonal matrix formed by extracting common row elements from G(S) 
so that A(S) contains numerator polynomial elements only. 


The decoupling pre-compensator is therefore given by: 


K(S) = A(S)" (eqn 3.2) 


However, this method is only possible if |A(s)| = O has all stable factors, 1e., ‘all 
roots of the charactoristic equation are in the left half plane. If this is not the case, 


then K(S) 1s formed from: 


Eqn 3.5 
Adjoint A(S) (eq 


K(S) = ae 
Realization Factor 
where the realization factor will contain the stable factors of |A(S)| together with other 
suitable lag elements to make K(S) physically realizable. 

The elements of the diagonal matrix R(S) and the elements of A(S) so formed are 
given in Table 4 and Table 5 respectively. As |A(S)| contains all stable factors, the 
ideal pre-compensator can be formed directly from Equation 3.2. The results of this 
Operation are given in Table 6. 

It is necessary to scale the columns of K(S) to ensure that the steady-state gains 
of the diagonal elements of Q(S) remain the same as those of G(S). The elements of 
this final pre-compensator are given in Table 7. Table 8 gives the speed related gain 
variation. 

Figure 3.2 and Figure 3.3 show that the addition of the pre-compensator has 


resulted in an improvement in outputs at 12 kts. 
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TABLE 4 


ELEMENTS OF DIAGONAL MATRIX R(S) 





i 
ame) = 2 
(Ss? + .06S + .25)(S* + .02S + .154)(8.2S + 1) 
1 
R,.(5) oon 
2 S(.434S + 1)(6.62S + 1)(4.18S + 1) 
| 
ate) = maT Avion 
(s* + .199S + .99)(10S + 1)(24S + 1) 
TABLE 5 
ELEMENTS OF MATRIX A(S) 
emis) = .25K,,(S° + .28 +.15)(8.25 + 1) 
ppGe)) = - 15K,,(S* + .06S + .25)(-8.57S + 1) 
Aso(S) = Kyo. 
yes) = -2K,,(105 + 1)(10.05s + 1) 
ese = K(8° + 22S + .1) 
emia = K,,(S* + .25°4 .1)(10705s + 1) 
byes) = A,,(S) = A,,(S) = 0 
A(S)| = .25K K..(S57 +.2S5 +.15)(S* + .28 + .1l 
Ee) | 11*22 333 23 + Ly(10. O58 + 1) : 





This “ideal” pre-compensator can be reduced to the individual elements of K(s) 
using step response data. 

After the reduction process is completed, the elements of the reduced pre- 
compensator so formed are given in Table 9. The terms K.,, K,,, and K,, given in 


Table 9 are the speed related compensator gain changes necessary for the pre- 
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TABLE 6 
ELEMENTS OF DIAGONALIZING PRE-COMPENSATOR 


_ . 00649 
~ 83.198? + 17.136S* + 9.58 + 1 
_ -.00454(-34.143s° + 28? - 8.295 +1) 
” 83. 19S? cssigenisiG. cee nie peel 
. 00285 
| .00376(10S + 1) 





C8) ~ “534° + 2788" + 182.6S° + 46.15S* + 11.58 + 1 
eee . 00013( 80050S° + 15650s* + 206658” + 2255S" + 3S + 1) 
ae 5352S° + 3326S° + 2113S* + 646S° + 16287 + 21.65 + 1 
P “a 00285 

ea 100S° +30S° +12S +1 


| Ki3(S) = K,,(S) = K3(S) = 0 


compensator to cope with the non-linearity of the warship model and these are the 


same as in Table 8. 
Figure 3.4 and Figure 3.5 show that interactions between inputs and outputs 


have not been eliminated in the system’s time response when we use the reduced order 


pre-compensator also. 
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TABRE 7 


FINAL ELEMENTS OF DIAGONALIZING PRE-COMPENSATOR 





i 

K S NS EEE 
1S) Bemose + 17.136S- + Sm5s + 1 

K.,(-34. 1438° + 2S? - 8.295 + 1) 
Kio(S) = 3 ES) 

Somes + 17. 136s” 2 o0Ss + 1 
K,,(S) = 1 

534S° + 2788S" + 182.6S° + 46.15S* + 11.5S + 1 
ms, = K.,(80050S° + 15650S* + 20665S° + 2255S? +35 + 1) 
af 5352S° + 3326S° + 2113S" + 646S° + 162S° + 21.6S + 1 
fens) = : 
33 H@Oc —ts0S- +126) cl 
Ki3(S) = Ky,(S) = K,,(S) = 0 


TABLE 8 
PRE-COMPENSATOR GAIN VARIATION 
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Figure 3.2 Step Response for Compensated Warship Model 
at Rudder Demand. 
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Figure 3.3 Step Response for Compensated Warship Model 
at Fin Demand. 
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FINAL REDUCED PRE-COMIVENSA TORE VE WIE. 





1 
Re 
5) 8.6S +l 
K.5(-34. 14S? + 25° - 8.35 + 1) 
RKi2( S) = Sey aD 
= §3.2S° “3-3 els S Veo Soe 1 
Re (sea 
RK..(S) = Bo2 
Os 14.8s* + 1.35 + 1) 
K ,(2000S* ="2s 471) 
RK,0(5) 0 = eee ee 
1308’? + 25S" + 12.65 + 1 
RK,4(S) = : 
$3 11.48 +1 
RK, ;(S) = RK,,(S) = RK,,(S) = 0 
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Figure 3.4 Step Response for Compensated Warship model 
at Rudder Demand using Reduced Order Pre-compensator. 
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Figure 3.5 Step Response for Compensated Warship model 
at Fin Demand using Reduced Order Pre-compensator. 
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IV. COMPUTER SIMULATION 


From the results discussed so far (all of which are from Roberts but were 
repeated at N.P.S. as part of this thesis), it can be seen that the pre-compensator 
substantially reduces the interaction intensity of the compensated system and results in 
improved ship stability and minimal loss of speed while executing normal maneuvers. 
Also, it has been shown that the outputs of ship’s roll angle, heading angle and speed 
can be controlled by a low-order pre-compensator. However, it can also be seen that 
large cross-coupling effects exist which cause long settling times. Therefore, a method 
must still be devised which controls the interactions between inputs and outputs and 
shortens settling time. The outputs of ship’s roll angle, heading angle and speed can be 
changed by a few factors. 

First of all, this thesis will introduce the philosophy of Function Minimization 
(F.M.), which is the main theme of this paper and introduce a few factors, which have 
important effects on outputs. 

Next, using above factors, this thesis will design other pre-compensators and 


simulate them. 


A. PHILOSOPHY OF FUNCTION MINIMIZATION 

Classical control theory has historically been applied to the design problem with 
the assistance of graphical presentations and trial and error methods. Such methods 
have been quite successful in the development of good control svstems, but do not 
answer the question “Is this the best system possible?”. 

If a given function has a minimum within the range of permitted parameter 
Variation, there exist numerical methods which can be used to find the minimum. 
These numerical methods have been programmed and most computer libraries contain 
one or more subroutines which can be used for Function Minimization. This thesis 
uses the HOOKE [Ref. 4] subroutine of DSL. 

Since one has the freedom to select the cost function to be used, certainly any of 
those used within optimal control theory can be chosen. We can therefore design an 
“optimal controller” without using the conventional theoretical approach. 

We can also use any of the well known performance indices as a cost function, 
1.e., NE at, fiE| dt, fiE|t dt are easily evaluated and minimized. 


oy 


In addition, one can select the cost function to suit the particular specification of 
the problem. 


Figure 4.1 shows the block diagram for simulation using Function Minimization. 
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Figure 4.1 Block Diagram of Function Minimization. 


B. CHOICE OF COMPENSATORS | 

Use of Function Minimization in the computer is relatively expensive, so some 
preliminary analysis and design is desirable to avoid excessive computer time. For a 
simple system ( one output ) a BODE design or a ROOT LOCUS design might be a 
good starting point. The compensator thus found is then simulated and improved 
through use of Function Minimization. For the ship problem in this thesis, the 
preliminary work was done by Roberts as can be seen in Chapters II and III. 
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This thesis uses Roberts’ compensators and outputs as a starting point. The 
desired output curves were obtained by modifving the known outputs. The questions 
to be answered are: 


e Can better performance be obtained using the original compensator but 
requiring a better output? 


e Can the compensator be reduced and the resulting output be as good or better 
than the original compensators? 


e Can a better output be obtained by increasing the order of the compensator? 

When the Function Minimization subroutine minimizes the function, it moves 
the poles and zeroes of the compensator to the best locations in the S-plane. If the 
compensator has too many poles and Zeroes, the program tends to set Z = P for the 
unneeded poles and Zeroes, thus the compensator order may be reduced. 

In many cases, even though the poles and zeroes may not exactly cancel, they are 
“placed so close to each other that thev contribute verv little to the result and so we 
may efficiently cancel them. Thus one can start the design by simply choosing each K. 
to have numerous Zeroes and poles as can be seen in Table 7 and Figure 3.2, which is 
the original pre-compensator. However, this thesis will show us that we don’t need 
that many poles and Zeroes for the choice of compensators. 

In this thesis, as shown in Table 9 and Figure 3.4, the reduced order pre- 
compensator was used. 

The optimum parameter values of this reduced order pre-compensator, which 
reduces the corresponding cross-coupling, can be determined by Function 
Minimization and Table 10 shows the optimum parameter values of the reduced order 
pre-compensator in Table 9 by Function Minimization. 

Figure 4.2 shows the positive damping effect on the system’s time response when 
optimum parameter values are used in the reduced pre-compensator. 

Figures 4.3 through Figure 4.5 show the difference between Roberts’ method and 
the Function Minimization method, where the continuous line is the output of Roberts 
and the dotted line is the output of Function Minimization. 

As shown in Figure 4.4 and Figure 4.5, cross-coupling is decreased dramatically 


in the case of Function Minimization. 


C. CHOICE OF DESIRED OUTPUTS 
One can always choose the ideal output as the desired output for the Function 


Minimization subroutine. However, this may not be a good choice because that 
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TABLE 10 


OPTIMUM PARAMETER VALUES OF REDUCED PRE- 
COMPENSATOR 


- 1 
$.65 qa 


K.o(-72. 148° + 1.98925S* - 2.38 + 1) 
76,20 ee? Coes ae le om ame 


7 . 63525 
32.0582 + 2.9758 # 1 
K.3(1741.258* = 1.943758 + 1) 
9.956" + Ti, 255° agro 
1 
See Ss 
= RK ( oye Peace 





particular output may be impossible to achieve. Although the minimization process 
will determine a closest fit, if the cost function is a “least squares” function the 
solution may not be acceptable. For example, with the ship problem of this thesis, the 
ideal output might be 

e Turn radius of two ship lengths 

e. Zero speed change 

e §«Zero roll angle 

None of these characteristics are possible. The results of a Function 

Minimization design will not satisfy any of them and the design achieved may not be 
acceptable. The desired output must therefore be chosen realistically, i.e., within the 
physical capabilities of the system. For the ship control problem, outputs were chosen 
based upon the results obtained by Roberts. In order to obtain better performance the 


desired outputs were chosen to be simular. 
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Figure 4.2 Output when we use the Optimum values of Reduced Pre-compensator 
at Rudder Demand. 
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Figure 4.3 Comparison between Roberts and F.M at Yaw output. 
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Figure 4.4 Comparison between Roberts and F.M at Roll output. 
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Figure 4.5 Comparison between Roberts and F.M at Speed output. 
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This thesis chose 3 desired system outputs , each consisting of a specification for 
heading, roll and speed. These desired output curves (AFGENI, 2 and 3) were 
obtained by modifying the known outputs (Figure 3.3) using Arbitrary Function 
Generators (AFGEN) provided by DSL. 

Figure 4.6, Figure 4.7 and Figure 4.8 show AFGENI, AFGEN2 and AFGEN3 
respectively. 

We can know that the AFGEN3 1s better than the other two cases by the trial 
and error method. Therefore this thesis will use the AFGEN3 for desired outputs for 


another simulation. 


D. COST FUNCTIONS 

When using Function Minimization as a design tool, a cost function must be 
chosen. This cost function is usually an integral. 

One possible procedure is to choose the desired performance as a reference and 
select a cost function which is the integral of the square of the difference between 
desired output and actual output. When the system has several outputs, the cost 


function must consider all of them, usually as a weighted sum. 


E. WEIGHTING FACTORS 
For the three outputs system, the cost function is of the form 
en 2 2 2 
Eee ie, dt + 1,/E* dt + A,)E,° dt. | 
There are no fixed rules for choosing the values factors A,, A, and 43. One 
approach is as follows: 


e Select the output which is considered most important, say # 1. Use this as a 
- reference and select A, = 1.0. 


e Base the value assigned to the second weighting factor, A, in this case, on the 
importance of output # 2 with respect to output # |. If output # 2 ts equally 
important as output # |, then A, = 1.0. If it is half as important, then 4, = 
=. 


e The third weighting factor is chosen in like manner. If the third output is less 
important than the first, then 4, < A, . If it is also less important than the 
second output, then A, < A,. 


The actual numbers chosen, for example h. = 0.1, are simply estimates based on 
experience. The designer may decide to change them after studying simulation results 
of a first design. 

In the case of ship control, the primary output is heading. Therefore, assign hi, a 
value of 1.0. Some change in speed is unavoidable and the desired response should 
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Figure 4.6 AFGEN lI. 
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Figure 4.7 AFGEN 2. 
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Figure 4.8 AFGEN 3. 
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Show this. However, if the desired response cannot be achieved perhaps more 
fluctuation in speed would be acceptable. For example, a value of approximately 0.1 
could be assigned to A.. If the desired roll output is twice as important as the speed 
output, then i) would be assigned a value of 0.2. 


Table 11 shows each case for simulation using weighting factors. 


TABLE 11 
SIMULATION CASES USING WEIGHTING FACTORS(W.F.) 


| | | 
CASE Le Ls | h, Sti. Urs 


L 1 mee bo | Figure 4.93 

. , 
2 il lou, Ge eed | Figure 4.10 | 
3 i —o 4 Figure 4.11 


aH 1 an i | Figure 4.12 


CASE4 is from Table 10, which uses Ay = |, h. = | and h. = | and this was 
repeated in Table 11 for comparison when we consider the weighting factors. 

Table 12 through Table 15 show the elements of the resulting compensators for 
each of the above when we used the F.M. approach. 

When F.M. was used for all coefficients in compensators, RK,,, RK,, and RRs 


were the parameters that the F.M. subroutine adjusted. 
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Figure 4.9 Step Response when A, = 1, Re)? and ae oet 


at Rudder Demand. 
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Figure 4.10 Step Response when A, = 1,4. =.6 and =.1 
at Rudder Demand. 
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Figure 4.11 Step Response when A, =1, A, =.6 and dh. =.4 
at Rudder Demand. 
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Figure 4.12 Step Response when A, = 1,4, =1 andd.=1 


at Rudder Demand. 
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TABEE 12 


OPTIMUM PARAMETEREVALE ES © ER er Dane. 
COMPENSAT ORVAT CASE 





1 
RK, ,($) = —————— 
Te 8.65 + 1 
K5(-47.64S° + 2.000625S° - 2.35 + 1) 
SE 
é 67.7S° 127. 38s + eos ar 
RiGee) eae 
gees 
RK, (3) = ———— 
20.4256 Bieec2ssue 1 
K.,(1915. 625s? - 2.001875S + 1) 
RR 50S) = ee eS ee 
90.625S° + 53.1258? + 23.855 + 1 
RK,,(S) = = 
a 11.45 + 1 


RK ;(S) = RK,,(S) = RK3(S) = 0 


ad 


TABLE 13 


Cee MOE RNIN Ener EL OES OF REDUCED PRE- 
COMPENSATOR AT CASE 2 W.F. 


1 
8.6S +l 
je 50 ot 1 99525c4g= 3235 + 1) 
45.25? + 31.13S* + 10.558 + 1 


. 603999 
25.8S* + 2.45 + 1 
Re@icoccemom 9755 + 1) 
53s? + g80s* + 34.65 + 1 
1 
11.48 + 1 
= RK,,(5) = RK;3(5) = 
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TABLE 14 


OPTIMUM PARAMETER VALUES OF REDUCED PRE- 
COMBIENSA [ORGAN H CA Silence.) aie 


1 
8.65 + 1 
Ks ( -5 6.067) + ores ez some) 
72S; +27. UsSe > eee meee 


6015 
25.357 + 2.355 +l 
K.,( 1842.58? - 1.9775S + 1) 
56.58? + 77.58? + 33.65 + 1 
_ 1 
eee ad 
= RK,,(S) = RK33(S) = 
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TABLE 15 
OPTIMUM PARAMETER VALUES OF REDUCED PRE- 


COMPENSATOR AT CASE 4 W.F. 





i 
eS 
t=) 8.6S +1 
K.o(-72.148° + 1.98925s* - 2.35 + 1) 
pe) oe Pe | 
POec oes 2SyOos + Vl SsS + 1 
Rieees = 1 
ee: = 63525 
al 32.058 + 2.9758 + l 
en aiez5o = 1, 943755 + 1) 
RK3.(S) = 3 a i ae 
9.258? + 111.258? + 47.15 + 1 


IL 
ivi oe eee 


= RK,.(S) = RK;,(S) = 0 





From simulation, it appears that the. most desirable output can be obtained at - 
ee A > .2andd>. = I. 

Figure 4.13 through Figure 4.15 show the comparison between the Figure 3.4 for 
Roberts, Figure 4.2 for CASE 4 and Figure 4.9 for CASE |. Where the output of 
Roberts is from Table 9, which gives the final reduced pre-compensator elements, the 
output of CASE 4 is from Table 10, which gives the optimum parameter values of 
reduced order pre-compensator by Function Minmization. The output of CASE | is 
the best case when we consider the weighting parameters for CASE I, 2 and 3 in Table 
fel 

The continuous line is used for Roberts, the dotted line for CASE 4 and the 
dashed line for CASE |! outputs. 

Figure 4.16 through Figure 4.18 show the comparison between the above three 
cases for Fin, Rudder and Power respectively. 
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Figure 4.13 Comparison between Fig 3.4, Fig 4.2 and Fig 4.9 
at Yaw Output. 
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Figure 4.14 Comparison between Fig 3.4, Fig 4.2 and Fig 4.9 
at Roll Output. 
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Figure 4.15 Comparison between Fig 3.4, Fig 4.2 and Fig 4.9 
at Speed Output. 
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Figure 4.16 Comparison between Roberts, CASE 4 and CASE | 
for Fin Output at 3° Rudder Demand. 
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Figure 4.17 Comparison between Roberts, CASE 4 and CASE | 
for Rudder Output at 3° Rudder Demand. 
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Figure 4.18 Comparison between Roberts, CASE 4 and CASE 1 
for Power Output at 3° Rudder Demand. 
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The simulation obtained so far was performed using three degrees for rudder 
demand and we can know that there is no problem about maximum values of input 
and output shown by Table 3 for each output. 

As can be seen in Figure 4.16 through Figure 4.18, CASE 1 is better than the 
other two cases. But we have to consider, carefully the physical realization in Figure 
4.18 because the power curve varies very fast. Therefore when we design the pre- 
compensator using weighting factors or other factors, we have to always consider 
physical realization limitations. 

Figure 4.19 through Figure 4.21 show the comparison between the above three 
cases for Fin, Rudder and Power when we use twenty degrees for rudder demand. 

From the above outputs, we can see that CASE 4 is better than the other two 


cases when we use twenty degrees for rudder demand. 


F. VARIATION OF NUMBER OF POLES AND ZEROES 

As shown by Table 7, the original pre-compensator by Roberts requires complex 
mathematical equations for a complete and detailed description. For many problems 
in the analysis and design of many modern dynamic systems, a simplified description, 
1.e., a low order model, is adequate and desirable. This thesis is concerned with the 
development of such low order models for the pre-compensator. | 

Two types of situations are commonly encountered in practice: 


e A system exists and can be tested, but its equations are not well known or not 
clearly defined. 


e A high order complex model of the system is Known and can be used, but 
it is undesirable for the problem to be studied. 


In either case the response of the system to a chosen signal can be obtained, and 
a low order model developed which has essentially identical outputs for identical inputs 
as the higher model. 

By carefully planned studies it should be possible to see how much reduction in 
order can be achieved, how closely the behavior of the.low order compares to that of 
the system and perhaps a best or optimum order can be found for the reduced order 
models. . 

In addition, by careful selection and classification of the pole-zero geometry of 
the high order system, it is hoped that a correlation may be found between such 
geometry and that of the best low order model. In any event, it is anticipated that 


some rules may be established for the choice of the number of poles and zeroes in the 
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Figure 4.19 Comparison between Roberts, CASE 4 and CASE | 
for Fin Output at 20° Rudder Demand. 


55 


60 


72 


48 36 64 


TIME (SEC) 


24 32 40 


16 


40 48 S56 64 72 80 
TIME (SEC) 


32 


FE 


~~ CITI ¥ISvI 


= 


~ 3389S) LISS 





Figure 4.20 Comparison betwecn Roberts, CASE 4 and CASE 1 
for Rudder Output at 20° Rudder Demand. 
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Figure 4.21 Comparison between Roberts, CASE 4 and CASE | 
for Power Output at 20° Rudder Demand. 
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low order model, such that evaluation of the pole, zero and gain will require only a few 
computer runs. 

As shown in Figure A.I of System Block Diagram for Simulation, K,, and K,, 
affect the Fin output of the system and K,,, K,, and K,, affect the Power output of 
the system. Therefore when we change the order of these elements, outputs, i.e., roll 
and speed will be changed by them also. 

Unfortunately there is no known mathematical basis for choosing the “best” 
order for a low order model for each element of the pre-compensator. For 
convenience, this thesis has changed the order of K,, and K,, for two cases only. 

Table 16 and Table 17 show the values of element for two cases and Figure 22 
through Figure 25 show the outputs for comparison between original reduced order 
pre-compensator by Table 10 and the above two cases respectively. Other elements are 


the same as in Table 10. 


TABI a6 


OPTIMUM PARAMETER VALUES OF REDUCED PRE- 
COMPENSATOR 
WHEN WE CHANGE K,, ONLY 


97.18759S + 3.94875 
scemiGGuigece + MOnScmeemIoOce 
K..(-72. 148° + 1.98925S* - 2.38 + 1) 
76.2S° + 28.63S* + 11.58 + 1 


. 63525 
Wesrtose- re eee 
K.3( 1741. 25S* - 1.943758 + 1) 
9.25S° + 111.25s* + 47.18 + 1 
i 
11.48 + 1 
Ko (5)! =" Xe ( Selo 
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TABLE 17 


OPTIMUM PARAMETER VALUES OF REDUCED PRE- 
COMPENSATOR 
WHEN WE CHANGE K,, AGAIN IN TABLE 16 


7 97.18759S + 3.94875 
7.53125S* + 21.075S + 12.025 
K.o(-72. 148° + 1.989258? - 2.35 + 1) 


76.2S° = 28.63S* + 11.5S + 1 


. 63525 
32.05S* + 2.9758 + l 
K.3( 1741. 258* - 1.94375S + 1) 
9.258? + 111.25S* + 47.18 + 1 
98.1875S + 2.05625 
6.65625S* + 28.05625S + 4.43625 
KS) = Kaas) = 0 





By the same algorithm, we can approach the appropriate order of pre- 
compensator and the values of each element for better design. 

As can be seen in Figure 4.22 through Figure 4.25, when we put the additional 

zero and pole to the element of the pre-compensator, the cross-coupling of the roll 


output is decreased dramatically. 
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Figure 4.22 Comparison between Table 10 and Table 16 
for Roll Output. 
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Figure 4.23 Comparison between Table 10 and Table 16 
for Speed Output. 
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Figure 4.24 Comparison between Table 10 and Table 17 
for Roll Output. 
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Figure 4.25 Comparison between Table 10 and Table 17 
for Speed Output. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In this thesis, the development of a pre-compensator to reduce the undesired and 
highly nonlinear cross-coupling effects in the maneuvering characteristics of a modern 
warship has been presented. Simulation results have shown that Function 
Minimization procedures for coordinated steering of a surface ship would significantly 
improve ship stabilitv, minimize loss of speed and reduce the interaction intensitv of 
the compensated system. 

Simulation results have also shown that the reduced pre-compensator can be 
determined by Function Minimization directly from an ideal pre-compensator based on 


given specifications. 


B. RECOMMENDATIONS 
Computer simulation for Function Minimization leads to the following 


recommendations: 


e This thesis uses the “HOOKE” subroutine of DSL for Function Minimization. 
It has certain constraint parameters, ie., ITMAX, CFT. Simulation outputs 
are changed by changing initial values and the step size for Function 
Minimization. Therefore, for given specifications initial values and step size for 
Function Minimization by constraint parameters should be determined by 
experience. 


e In this thesis, it has been shown that Function Minimization can be used for 
various cases as shown in Table 11. A particular case for needed specification 
should be determined and simulated by trial and error method for element 
optimum vaiues which have undesired cross-coupling effects. 


e In this thesis, it has been assumed that the stabilizer fin to yaw cross-coupling 
term, G,,(S), is a null entry for the class of warship considered. However, it is 
recommended that cross-coupling between these parameters be considered in 
further studies on general surface ships. 


e Further research should investigate ship characteristic constraints. 


e Further research should investigate the effects of various sea state conditions, 
speeds and maneuvering. 
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APPENDIX A 
SYSTEM BLOCK DIAGRAM FOR SIMULATION 


RIN 


SIN 





Figure A.1 System Block Diagram for Simulation. 
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APPENDIX B 
CONSTRAINT PARAMETERS OF FUNCTION MINIMIZATION 


To use HOOKE, the user must initialize the following arguments and, in the 
main program, 
CALL HOOKE (X, STEP, N, ITMAX, CFYOL, ALPHA, BETA, Cr aga ee: 
W, IPRINT, MINMAX) 
All of these arguments must be initialized in MAIN, except for X, CF, Q, QQ, 
and W. Recommended values for ALPHA and BETA are ALPHA = 2., BETA = 0.5. 


All of the arrays, 1.e., X, STEP, Q, QQ and W, must be declared and dimensioned in 
MAIN. 
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TABLE 18 


THE LIST OF PARAMETERS 





| ARGUMENT MEANING 


the aera y of N parameter values. The user must 
e 


x Supply. & 


initial guesses, either in the DSL 


program or in MAIN. 


an array of dimension N containin a Pad et ae 
| See stepsizes to be used in the searc 





at most 


the Bue rues parameters (a positive integer, 


ithe maximum number of function calls to be 


ITMAX performed. 


the error in the criterion function to be 

ee TOL reached pefore the program terminates 
(difference between the current value and the 
Be cVa One: stage value). 


the factor of (Y =- X) which is added to a to get 
ALPHA ANEW; a number greater than or equal to 


oo ee reduction factor; a number between 
_— 


‘cr the value of the criterion function. 


: 





arrays of dimension N, to be used as work Bec e: 
2700, They must be declared and dimensioned in t! 


program. 


an integer fla 
IPRINT 2 a 


QO for no intermediate printout 
ior intermediate printout 


Of xX, CE, the number of function evaluation and 
step-reduction. 


notification of 


an integer flag: 
MINMAX 


-l1 searches for a minimum 
+1 searches for a maximum 


67 





APPENDIX C 
COMPUTER PROGRAM FOR UNCOMPENSATED SYSTEM 


TITLE SIMULATION OF UNCOMPENSATED SYSTEM 
ARRAY Al(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4(4).... 
A5(1),B5(3),46(1),B6(2) 
TABLE Al(1)= 1,B1(1-3) = 4,.24,1,A2(1-2) = -8.57,1,B2(1-4) = 53.3,17.17.... 
9,52,1,A3(1)= 1,B3(1-5)= 12,32.25,11.2,1,0,A4(1-2)= 10, 1.... 
B4( 1-4) = 240,58,26,1,45(1) = 1,B5(1-3) = 240,34,1,A6(1) = 1... 
B6(1-2)= 24,1 
CONST K11=0.114,K12=0.18,K22 = .01,K31 = 0.058,K32 = 0.096,K33 = 0.1... 
RL=0.,YW=0.,SP=0. 
DERIVATIVE 
FIN =0.0*STEP(0) 
RUDDER = 3.*STEP(0) 
POWER =0.*STEP(0) 
ROLLI = TRNFR(0,2.RL,A1.B1,K11*FIN) 
ROLL2= TRNFR(1,3,RL,A2.B2.K12*RUDDER) 
ROLL = ROLL1+ ROLL2 
YAW =TRNFR(0,4, YW,A3,B3,K22*RUDDER) 
SPEED1 = TRNFR(1,3,SP,A4,B4,K 31*FIN) 
SPEED2= TRNFR(0,2,SP,A5,B5,K32*RUDDER) 
SPEED3= TRNFR(0,1,SP,A6,B6,K33*POWER) 
SPEED = -SPEEDI-SPEED2 + SPEED3 
CONTROL FINTIM = 80. 
*RINT 1.,ROLL,YAW,SPEED 
SAVE (SI) 0.1,ROLL,YAW,SPEED 
GRAPH(G1/S1,DE= TEK618,PO=O0,.5) TIME(LE=8.,UN= SEC) ROLL(LO=-1,LI=.... 
SC = .25,NI=8), YAW(LO=-1,LI = 3,SC =.25,NI=8).... 
SPEED(LO=-1,LI=4,SC =.25,NI = 8) 
LABEL(G1)STEP RESPONSE FOR UNCOMPENSATED WARSHIP MODEL AT 12 KTS 
LABEL(G1)10 % RUDDER DEMAND 
END 
STOP 
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APPENDIX D 
COMPUTER PROGRAM FOR ORIGINAL PRE-COMPENSATOR 


TITLE SIMULATION OF ORIGINAL PRE-COMPENSATOR 
ARRAY AIl(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4(4).... 
A5(1),B5(3),A6(1),B6(2) 
ARRAY C1(1),D1(4),C2(4),D2(4),C3(2),D3(6),C4(6),D4(7),C5(1),D5(4) 
TABLE Al(1)=1,B1(1-3)= 4,.24, 1,4 2(1-2) = -8.57,1,B2(1-4) = 53.3,17.17,... 
9.52,1,43(1) = 1,B3(1-5) = 12,32.25,11.2,1,0,44(1-2) = 10,1... 
B4( 1-4) = 240,58,26,1,A.5(1) = 1,B5(1-3) = 240,34,1,A6(1)= 1... 
B6(1-2) = 24.1 
TABLE Cl(1)=1,D1(1-4) = 53.19,17.136,9.5,1,... 
C2(1-4)= -34.143,2,-8.29, I... 
D2(1-4) = 53.19,17.13,9.5,1,C3(1-2) = 10,1... 
D3(1-6) = 534,278, 182.6,46.15,11.5,1,... 
C4(1-6) = 80050, 15650,20665,2255,3,],... 
Dd( 1-7) = 5352,3326.21 13,646, 162,21.6,1,... 
C5(1)= 1,D5(1-4) = 100,30, 12,1 
CONST K11=0.114,K12=0.18,K22 =.01,K31 =0.058,K32=0.096,K33=0.1,... 
RL=0.,YW=0.,SP=0. | 
CONST GC10=0,GC20 =0,GC30=0,GC40=0,GC50=0 
DERIVATIVE | 
RIN =0.0*STEP(0) 
YIN=3.*STEP(0) 
SIN =0.*STEP(0) 
El=RIN-ROLL 
E2=YIN-YAW 
E3=SIN-SPEED 
FINI = TRNFR(0,3,GC10,C1,D1,E1) 
FIN2= TRNFR(3,3,GC20,C2,D2,-1.6*E2) 
FIN= FINI + FIN2 
ROLLI=TRNFR(0,2,RL,A1,B1,K11*FIN) 
RUDDER = YIN-YAW 
ROLL2= TRNFR(1,3,RL,A2,B2,K12*RUDDER) 
ROLL=ROLL1+ROLL2 
YAW =TRNFR(0,4, YW,A3,B3,K22*RUDDER) 
SP1 = TRNFR(1,5,GC30,C3,D3,0.579*E1) 
SP2= TRNFR(5,6,GC40,C4,D4,0.0456*E2) 
SP3= TRNFR(0,3,GC50,C5,D5,E3) 
SPEED 1 = TRNFR(1,3,SP,A4,B4,K31*FIN) 
SPEED2= TRNFR(0,2,SP,A5,B5,K32*RUDDER) 
POWER = SP3-SP2-SP1 
SPEED3 = TRNFR(0,1,SP,A6,B6,K33*POWER) 
SPEED = SPEED3-SPEED2-SPEED1 
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CONTROL FINTIM = 80. 
*RINT L.,ROLL,YAW,SPEED 
SAVE (S1) 0.1,ROLL,YAW,SPEED 
GRAPH(G1/S1,DE= TEK618,PO=0,.5) TIME(LE=8.,UN = SEC) ROLL(LO=-2,LI =... 
SC =.5,NI=8), YAW(LO=-2,LI = 3,SC=.5,NI = 8)... 
SPEED(LO = -2,LI = 4,SC =.5,NI = 8) 
LABEL(G1)STEP RESPONSE FOR COMPENSATED WARSHIP MODEL AT 12 KTS 
LABEL(G1)10 % RUDDER DEMAND 
END 
STOP 
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APPENDIX E 


COMPUTER PROGRAM FOR ORIGINAL REDUCED PRE- 
COMPENSATOR 


TITLE SIMULATION FOR ORIGINAL REDUCED ORDER PRE-COMPENSATOR 
ARRAY AI(1),B1(3),A2(2),B2(4),A 3(1),B3(5),A4(2),B4(4),... 
A35(1),B5(3),A6(1),B6(2) 
ARRAY C1(1),D1(2),C2(4),D2(4),C3(1),D3(3),C4(3),D4(4),C5(1),D5(2) 
TABLE Al(1)= 1,B1(1-3) = 4,.24, 1, A2( 1-2) = -8.57,1,B2( 1-4) = 53.3,17.17.... 
9.52,1,A3(1) = 1,B3(1-5) = 12,32.25,11.2,1,0,A4(1-2) = 10,1... 
Ba( 1-4) = 240,58,26,1,A.5(1) = 1,B5(1-3)= 240,34,1,A6(1) = 1,... 
B6(1-2)= 24,1 
TABLE Cl(1)= 1,D1(1-2) = 8.6,1,C2(1-4) = -34.14,2,-3.3,1,... 
D2(1-4) = 53.2,17.13,9.5,1,C3(1) = 1,D3(1-3) = 14.8,1.3,1,... 
C4(1-3) = 2000,-2,1,D4(1-4) = 130.25,12.6,1,... 
C5(1)= 1,D5(1-2) = 11.4,1 
CONST K11=0.114,K12=0.18,K22=.01,K31=0.058,K32=0.096,K33 =0.1.... 
RL=0.,YW=0.,SP=0. 
CONST KC1=-1.6,KC2=0.579,KC3 = 0.0456... 
GC10=0,GC20= 0,GC30 = 0,GC40 = 0,GC50=0 
DERIVATIVE 
RIN =0.0*STEP(0) 
YIN = 3.*STEP(0) 
SIN = 0.*STEP(0) 
El=RIN-ROLL 
E2= YIN-YAW 
E3 = SIN-SPEED 
FINI =TRNFR(0,1,GC10,C1,D1,E1) 
FIN2= TRNFR(3,3,GC20,C2,D2,KC1*E2) 
FIN= FINI + FIN2 
ROLL! = TRNFR(0,2,RL,A1,BI,K11*FIN) 
RUDDER= YIN-YAW 
ROLL2= TRNFR(1,3,RL,A2,B2,K12*RUDDER) 
ROLL =ROLL1I+ROLL2 
YAW =TRNFR(0,4,YW,A3,B3,K22*RUDDER) 
SP1 = TRNFR(0,2,GC30,C3,D3,KC2*E1) 
SP2 = TRNFR(2,3,GC40,C4,D4,KC3*E2) 
SP3 = TRNFR(0,1,GC50,C5,D5,E3) 
SPEED1 = TRNFR(1,3,SP,A4,B4,K31*FIN) 
SPEED2 = TRNFR(0,2,SP,A5,B5,K32*RUDDER) 
POWER = SP3-SP2-SP1 
SPEED3 = TRNFR(0,1,SP,A6,B6,K33*POWER) 
SPEED = SPEED3-SPEED2-SPEED! 
CONTROL FINTIM = 80. 
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*RINT 1.,ROLL,YAW,SPEED 
*RINT 1.,FIN,RUDDER,POWER 
SAVE (S1) 0.1,ROLL,YAW,SPEED 
GRAPH(G1/S1,DE = TEK618,PO = 0,.5) TIME(LE=8.,UN= SEC) ROLL(LO=-2,LI=1.... 
$C =.5,NI= 8), YAW(LO=-2,LI =3,SC=.5,NI= 8)... 
SPEED(LO = -2,LI =4,SC =.5,NI = 8) 
LABEL(G1)STEP RESPONSE FOR COMPENSATED WARSHIP MODEL AT 12 KTS 
LABEL(GI)FIGURE 3.4 
END 
STOP 


IZ 


APPENDIX F 
COMPUTER PROGRAM FOR AFGEN SUBROUTINE 


TITLE SIMULATION FOR desired outputs using afgen subroutine 
AFGEN RL=0,0,1,-.03,4,-.1,10,0,16,0,23,0,29,0,36,0.... 
42,0,48,0,54,0,61,0,67,0,80,0 
AFGEN YW =0,0,1,.00006,4,.006,8,.04,10,.07,20,.29,30,.55,40,.8,80,1.61 
AFGEN SD=0,0,1,-.004,10,-.01,21,0,31,0,41,0,51,0,70,0,80,0 
DERIVATIVE 
X=TIME 
ROLL = NLFGEN(RL,X) 
YAW = NLFGEN(YW,X) 
SPEED = NLFGEN(SD,X) 
CONTROL FINTIM = 80 
Sein] 1.,ROLL,YAW,SPEED 
SAVE (S1) 0.1.ROLL,YAW,SPEED | 
GRAPH(G1/S1,DE= TEK618,PO=0,.5) TIME(LE=8.,UN = SEC), ROLL(LO=-?2,LI = 1.... 
SC =.5,NI=3), YAW(LO= -2,LI =3,SC=.5,NI = 8)... 
SPEED(LO= -2,LI=4,SC=.5,NI= 8) 
LABEL(G1) AFGN 1 FOR FIGURE 4.6 
END 
STOP 
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APPENDIX G 
COMPUTER PROGRAM FOR FUNCTION MINIMIZATION 


D COMMON/HANDSJ/FLAG,ER,K1,K2,K3,K4,K5,K6,K7,K8,K9,K10,K11,K12 | 
D COMMON/HANDJ/K13,K14 | 
TITLE SIMULATION FOR FUNCTION MINIMIZATION SUBROUTINE } 
ARRAY A1(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4(4),A5(1),B5(3),... | 
A6(1),B6(2) 
STORAG C1(1),D1(2),C2(4),D2(4),C3(1),D3(3),C4(3),D4(4),C5(1),D5(2) 
TABLE Al(1)=1,B1(1-3)= 4,.24, 1,A2(1-2) = -8.57,1,B2(1-4) = 53.3,17.17,... | 
9.52,1,A3(1)= 1,B3(1-5)= 12,32.25,11.2,1,0,A4(122}= low... | 
B4(1-4) = 240,58,26,1,A5(1) = 1,B5(1-3) = 240,34, 1... | 
A6(1)= 1,B6(1-2) = 24,1 
CONST K01=-34.14.K20= 2,K30= -8.3,K40= 53.2,K50= 17.13,K60 = 9.5,K70=.579.... | 
K80= 14.8.K90= 1.3,K 100 = 2000,K110 = -2,K120 = 130,K 130 = 25... | 
K140= 12.6 , 
CONST Cll =.114.C12=.18,C22 = .01,C31 = .058,C32 = .096,C33 = .1,IC =0.,... 
KC1=-1.6,KC3 = 0.0456 
PARAM KIMIN=-50,KIMAX = 200, K2MIN = .01,K2MAX = 200, K3MIN =-10,K3MAX = 200,, 
K4MIN=.01,K4MAX = 200, KSMIN=.01,K5MAX = 200, K6MIN=.01,K6MAX = 200... 
K7MIN=.01,K7MAX = 200, K8MIN =.01,K8MAX = 200, K9MIN = .01,K9MAX = 200... 
KIOMIN=.01,KIOMAX = 3000,K11IMIN=-10,KI1IMAX = 200,K12MIN=.01.... 
K12MAX = 200,K13MIN=.01,KI3MAX = 200,K14MIN =.01,K14MAX = 200 
AFGEN RL=0,0,1,-.03,4,-.1,10,0,16,0,23,0,29,0,36,0,... . 
42,0,48,0,54,0,61,0,67,0,80,0 
AFGEN YW=0,0,1,.00006,4,.006,8,.04, 10,.07,20,.29,30,.55,40,.8,80,1.61 
AFGEN SD =0,0,1,-.004,10,-.01,21,0,31,0,41,0,51,0,70,0,80,0 
INITIAL SEGMENT 
IF(FLAG.LT.0O.) K1 = KOl 
IF(FLAG.LT.0.) K2= K20 
[F(FLAG.LT.0.) K3 = K30 
[F(FLAG.LT.0.) K4= K40 
[F(FLAG.LT.0.) K5=K50 
[F(FLAG.LT.0.) K6= K60 
IF(FLAG.LT.0.) K7= K70 
IF(FLAG.LT.0.) K8= K80 
[F(FLAG.LT.0.) K9= K90 
IF(FLAG.LT.0.) K10= K100 
IF(FLAG.LT.O.) K11=K110 
IF(FLAG.LT.0.) K12=K120 
I[F(FLAG.LT.0O.) K13= K130 
IF(FLAG.LT.O.) K14=K140 
FLAG=FLAG#1 
IF((K1.LE.KIMIN).OR.(K1.GE.KIMAX)) THEN 


. 
| 
! 
| 
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K1=KOl 

ENDIF 

IF((K2.LE.K2MIN).OR.(K2.GE.K2MAX)) THEN 
K2=K20 

ENDIF 

IF((K3.LE.K3MIN).OR.(K3.GE.K3MAX)) THEN 
K3= K30 

ENDIF 

IF((K4.LE.K4MIN).OR.(K4.GE.K4MAX)) THEN 
K4=K40 

ENDIF 

[F((K5.LE.KSMIN).OR.(K5.GE.KSMAX)) THEN 
K5=KS0 

ENDIF 

[F((K6.LE.K6MIN).OR.(K6.GE.K6MAX)) THEN 
K6=K60 

ENDIF 

[F((K7.LE.K7MIN).OR.(K7.GE.K7MAX)) THEN 
K7=K70 

ENDIF 

[F((K8.LE.K8MIN).OR.(K8.GE.K8MAX)) THEN 
K8 = K80 

ENDIF 

IF((K9.LE.K9MIN).OR.(K9.GE.K9IMAX)) THEN 
K9 = K90 

ENDIF | 

IF((K10.LE.K1OMIN).OR.(K10.GE.K10MAX)) THEN 
K10= K100 

ENDIF | 

[F((K11.LE.K11MIN).OR(K11.GE.K11MAX)) THEN 
KI1=K110 

ENDIF 

IF((K12.LE.K12MIN).OR.(K12.GE.K12MAX)) THEN 
K12=K120 

ENDIF 

IF((K13.LE.KI3MIN).OR.(K13.GE.K13MAX)) THEN 
K13= K130 

ENDIF 

[F((K14.LE.K14MIN).OR.(K14.GE.K 14MAX)) THEN 
K14=K140 

ENDIF 

C1(1)=1 

DI(1)=8.6 

D1(2)=1 

C2(1)=K1 

C2(2)= K2 

C2(3)=K3 
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C2(4)= 1 
D2(1)=K4 
D2(2)= KS 
D2(3)= K6 
D2(4)= 1 
C3(1)= K7 
D3(1)=K8 
D3(2)=K9 
D3(3)=1 
C4(1) = K10 
C4(2)=K11 
C4(3)= 1 
D4(1)= K12 
D4(2) = K13 
D4(3)=K14 
D4(4)=1 
C5(1)=1 
DS5(1)=11.4 
DS(2)=1 
DERIVATIVE 
RIN =0.0*STEP(0) 
YIN = 3.*STEP(0) 
SIN =0.*STEP(0) 
El=RIN-ROLL 
E2= YIN-YAW | 
E3=SIN-SPEED 
FINI =TRNFR(0,1,IC,C1,D1,E1 
FIN2 = TRNFR(3,3,IC,C2,D2,KC1*E2) 
FIN =FINI+ FIN2 
ROLL1 = TRNFR(0,2,1C,A1,B1,C11*FIN) 
RUDDER= YIN-YAW 
ROLL2= TRNFR(1,3.1C,A2,B2,C12*RUDDER) 
ROLL = ROLL1+ ROLL? 
YAW = TRNFR(0,4,IC,A3,B3,C22*R UDDER) 
SP1=TRNFR(0,2,IC,C3,D3,E1) 
SP2= TRNFR(2,3,IC,C4,D4,KC3*E2) 
SP3 = TRNFR(0,1,1C,C5,DS5,E3) 
SPEED1 = TRNFR(1,3,1C,A4,B4,C31*FIN) 
SPEED2 = TRNFR(0,2,IC,A5,B5,C32*RUDDER) 
POWER = SP3-SP2-SP1 
SPEED3 = TRNFR(0,1,1C,A6,B6,C33*POWER) 
SPEED = SPEED3-SPEED2-SPEED1 
X=TIME 
YNL1 = NLFGEN(RL,X) 
YNL2= NLFGEN(YW,X) 
YNL3 = NLFGEN(SD,X) 
EE1 =((ROLL-YNL1)**2) 
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El11=INTGRL(IC,EE1) 
EE2=((YAW-YNL2)**2) 
E22= INTGRL(IC,EE2) 
EE3 = ((SPEED-YNL3)**2) 
E33 = INTGRL(IC,EE3) 
E=E11+E22+E33 


CONTROL FINTIM = 50,DELT=.01 
TERMINAL 


eae 


FORTRAN 


[IMPLICIT REAL*8(A-H,O-Z) 
DIMENSION X(14),STEP(14),Q(14),QQ(14),W(14) 
STEP(1)=4. 

STEP(2)=.001 

STEP(3) =2. 

STEP(4) = 4. 

STEP(5) = 2. 

STEP(6) = 1. 

STEP(7) = 0.005 

STEP(8)= 1. 

STEP(9)=0.1 

STEP(10)= 15. 

STEP(11) = .005 

STEP(12)=7. 

STEP(13)=5. 

STEP(14) =2. 

N=14 

ITMAX = 200 

CFTOL = .0000001 

ALPHA =2 

BETA =.5 

IPRINT=0 

MINMAX =-1 

CALL HOOKE(X,STEP,N,ITMAX,CFTOL,ALPHA,BETA, 
*CF,Q,Q0Q,W,IPRINT,MINMAX) 
STOP 


~ END 
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APPENDIX H 


COMPUTER PROGRAM FOR COMPARISON BETWEEN ROBERTS 
AND F.M. 


TITLE SIMULATION FOR COMPARISON BETWEEN ROBERTS AND F.M. 
*#*4% ROBERTS OUTPUT ****# 
ARRAY AI1(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4(4).... 
A5(1),B5(3),A6(1),B6(2) 
ARRAY C1(1),D1(2),C2(4),D2(4),C3(1),D3(3),C4(3),D4(4),C5(1),D5(2) 
TABLE Al(1)=1,B1(1-3) = 4,.24,1,A2(1-2) = -8.57,1,B2(1-4) = 53.3,17.17,... 
9.52,1,A3(1) = 1,B3(1-5) = 12,32.25,11.2,1,0,A4(1-2) = 10,1,... 
B4(1-4) = 240,58,26,1,4.5(1) = 1,B5(1-3) = 240,34,1,A6(1) = 1... 
B6(1-2) = 24,1 
TABLE C1(1)=1,D1(1-2)= 8.6,1,C2(1-4) = -34.14,2,-8.3,1,... 
D2( 1-4) = 53.2,17.13,9.5,1,C3(1) = 1,D3(1-3) = 14.8, 1.3,1,... 
C4(1-3) = 2000,-2,1,D4(1-4) = 130,25, 12.6,1,... 
C5(1)= 1,D5(1-2) = 11.41 
CONST K11=0.114,K12=0.18,K22 = .01,K31 = 0.058,K32 = 0.096,K33 = 0.1... 
RL=0.,YW=0.,SP=0. _ 
CONST KC1=-1.6,KC2=0.579,KC3 = 0.0456... 
GC10=0,GC20=0,GC30 = 0,GC40 = 0,GC50=0 
DERIVATIVE 
RIN = 0.0*STEP(0) 
YIN = 3.*STEP(0) 
SIN = 0.*STEP(0) 
El=RIN-ROLL 
E2=YIN-YAW 
E3=SIN-SPEED 
FINI = TRNFR(0,1,GC10,C1,D1,E1) 
FIN2= TRNFR(3,3,GC20,C2,D2,K C1*E2) 
FIN = FINI + FIN2 
ROLL1=TRNFR(0,2,RL,A1,B1,K11*FIN) 
RUDDER= YIN-YAW 
ROLL2=TRNFR(1,3,RL,A2,B2,K 12*RUDDER) 
ROLL=ROLL1+ROLL2 
YAW = TRNFR(0,4,YW,A3,B3,K22*RUDDER) 
SP1 = TRNFR(0,2,GC30,C3,D3,KC2*E1) 
SP2 = TRNFR(2,3,GC40,C4,D4,KC3*E2) 
SP3 = TRNFR(0,1,GC50,C5,D5,E3) 
SPEED1 = TRNFR(1,3,SP,A4,B4,K31*FIN) 
SPEED2= TRNFR(0,2,SP,A5,B5,K32*RUDDER) 
POWER = SP3-SP2-SP1 
SPEED3 = TRNFR(0,1,SP,A6,B6,K33*POWER) 
SPEED = SPEED3-SPEED2-SPEEDI 


78 


**%%** ELINCTION MINIMIZATION OUTPUT ***** 
ARRAY AAI(1),BB1(3),AA2(2),BB2(4),AA 3(1),BB3(5),AA4(2),BB4(4).... 
AAS(1),BB5(3),AA6(1),BB6(2) 
ARRAY CC1(1),DD1(2),CC2(4),DD2(4),CC3(1),DD3(3) 
ARRAY CC4(3),DD4(4),CC5(1),DD5(2) 
TABLE AA1(1)= 1,BB1(1-3)= 4,.24,1,AA2(1-2) = -8.57, 1... 
BB2(1-4) = 53.3,17.17,9.52,1,... 
AA3(1) = 1,BB3(1-5) = 12,32.25,11.2,1,0,AA4(1-2) = 10,1... 
BB4(1-4) = 240,58,26,1,AA5(1)= 1,BB5(1-3) = 240,34,1,AA6(1)= 1... 
BB6(1-2)= 24,1 
TABLE CCI(1)= 1,DD1(1-2) = 8.6,1,CC2(1-4) = -72.14,1.98925,-2.3, 1... 
DD2(1-4) = 76.2,28.63,11.5.1,DD3(1-3) = 32.05,2.975. 1... 
CC4(1-3) = 1741.25,-1.94375,1,DD4(1-4) = 9.25,111.25,47. 1,1... 
CC5(1)= 1,.DD5(1-2) = 11.4,1,CC3(1) = .63525 
CONST KK11=0.114,KK12 =0.18,KK22=.01,KK31 =0.058,KK32=0.096,KK33=0.1.... 
RRL=0.,YYW =0.,SSP=0.,KKC1 =-1.6, KKC3 = 0.0456 
CONST GGC10=0,GGC20 = 0,GGC30= 0,GGC40 = 0,GGC50=0 
DERIVATIVE 
RRIN =0.0*STEP(0) 
YYIN = 3.*STEP(0) 
SSIN =0.*STEP(0) 
EE1 = RRIN-RROLL 
EE2= YYIN-YYAW 
EE3=SSIN-SSPEED 
FFIN1 = TRNFR(0,1,GGC10,CC1,DD1,EE1) 
FFIN2 = TRNFR(3,3,GGC20,CC2,DD2,KKC1*EE2) 
FFIN = FFIN1 + FFIN2 
RROLLI=TRNFR(0,2,RRL,AA1,BB1,KK11*FFIN) 
RRUDDE= YYIN-YYAW 
RROLL2=TRNFR(1,3,RRL,AA2,BB2,K K12*RRUDDE) 
RROLL=RROLL1+RROLL?2 
YYAW = TRNFR(0,4, YYW.AA3,BB3,KK22*RRUDDE) 
SSP1 = TRNFR(0,2.GGC30,CC3,DD3,EE1) 
SSP2 = TRNFR(2,3,GGC40,CC4,DD4,.KK C3*EE2) 
SSP3 = TRNFR(0,1,GGC50,CC5,DDS5,EE3) 
SSPEE1 = TRNFR(1,3,SSP,AA4,BB4,K K31*FFIN) 
SSPEE2 = TRNFR(0,2,SSP,AA5,BB5,KK32*RRUDDE) 
PPOWER = SSP3-SSP2-SSP1 
SSPEE3 = TRNFR(0,1,SSP,AA6,BB6,KK33*PPO WER) 
SSPEED = SSPEE3-SSPEE2-SSPEE1 
CONTROL FINTIM = 80. 
*RINT 1.,FFIN,RRUDDER,RROLL,YYAW,SSPEED 
SAVE (S1) 0.1,ROLL,RROLL,YAW,YYAW,SPEED,SSPEED 
GRAPH(G1/S1,DE = TEK618,PO=0,.5) TIME(LE = 8.,UN = SEC) YAW(LO=-2,LI=1.... 
SC =.5,NI=8), YYAW(LO= -2,LI1 =4,SC = .5,NI = 8) 
LABEL(G1)COMPARISON BETWEEN ROBERT AND F.M 
LABEL(G1)FOR YAW OUTPUT 
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GRAPH(G2/S1,DE = TEK618,PO=0,.5) TIME(LE=8.,UN = SEC) ROLL(LO=-2,LI=1.... 
SC=.5,NI=8),RROLL(LO=-2,LI=4,SC=.5,NI=8) 

LABEL(G2)COMPARISON BETWEEN ROBERT AND F.M 

LABEL(G2)FOR ROLL OUTPUT 

GRAPH(G3/S1,DE = TEK618,PO=0,.5) TIME(LE=8.,UN = SEC) SPEED(LO = -2,LI= 1... 
SC =.5,NI = 8),SSPEED(LO = -2,LI = 4,SC = .5,NI = 8) 

LABEL(G3)COMPARISON BETWEEN ROBERT AND F.M 

LABEL(G3)FOR SPEED OUTPUT 

END 

STOP 
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